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INTRODUC  TI  ON 


A  program  of  applied  research  in  the  problem  of  aerodynamic 
heating  of  guided  missiles  has  been  in  progress  as  a  funded  task  at  the 
Naval  Ordnance  Test  Station  for  the  past  three  years.  In  addition  to 
several  relatively  long-term  major  efforts,  a  number  of  investigations 
of  smaller  scope  form  a  significant  portion  of  this  prograpi.  In  many 
instances,  these  latter  investigations  grew  out  of  or  were  extensions 
to  work  required  for  a  specific  development  program.  In  these  cases, 
the  results  were  generally  published  in  local  informal  reports  of  very 
limited  distribution,  or  were  submerged  in  a  more  extensive  report 
concerning  the  related  development  task.  • 

• 

This  report  is  one  of  sev'eral  that  are  being  published  to  provide 
easier  access  and  wider  distribution  to  the  results  attained  in  these 
investigations.  These  reports  consist  of  a  collection  of  notes  on  related 
topics  gleaned  from  the  informal  reports  which  were  published  for  local 
distribution  only. 

The  notes  in  this  report  are  concerned  with  analytical  techniques. 
The  methods  described  utilize  analog  computers,  and  criteria  are  given 
under  which  more  simplified  methods  may  be  profitably  employed. 
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SECTION  I.  TEMPERATURE  CHANGES  DURING  HEATING 
OF  AN  HPAG  5.  0  INCH  MOTOR: 
CORRELATION  OF  EXPERIMENTAL  DATA 
WITH  ANALOG  COMPUTATIONS 

By 

C.  H.  Johnson 


INTRODUCTION 

A  major  difficulty  in  estimating  temperature  histories  inside  missile 
components  subjected  to  unsteady-state  aerodynamic  heating  is  in 
setting  up  a  suitable  model  and  finding  accurate  physical  property  data 
to  use  in  the  computations.  Deficiencies  in  thermal  property  data  for 
engineering  materials,  not  only  with  respect  to  accuracy  at  ambient 
temperatures  but  also  in  regard  to  temperature  dependence,  make  it 
impractical  to  consider  any  other  than  principal  effects  on  heat  transfer 
in  setting  up  the  analytical  model.  Nevertheless,  it  is  important  for 
good  simulation  to  include  all  the  essential  effects,  and  laboratory  ex¬ 
periments  may  be  needed  to  determine  which  effects  should  be  considered 
in  analytical  studies.  This  note  illustrates  how  temperature  data  on  a 
5.  0-inch  HPAG  rocket  motor  served  to  guide  the  set-up  of  the  analog 
computer  model.  The  study  was  prompted  by  observations  that  temper¬ 
atures  in  Sidewinder  1C  motors,  heated  in  the  NOTS  radiant  heat 
facility,  differed  appreciably  depending  on  whether  the  test  point  was 
at  the  top,  side,  or  bottom  of  the  motor. 


THE  PROBLEM 

The  HPAG  motor  was  used  as  the  experimental  model  in  this  study 
because  of  its  similarities  to  current  missile  propulsion  units.  Data 
on  this  and  other  motors  subjected  ‘o  oven  heating  and  cooling  tests  are 
given  in  Ref.  1.  The  present  study  is  based  on  a  step  increase  in 
temperature,  78®  to  140 *F.  This  interval  was  chosen  so  that  available 
thermal  conductivity,  specific  heat,  and  thermal  expansion  data  would 
apply,  and  the  effect  of  temperature  on  these  properties  would  be 
minimized.  All  motor  grain  physical  properties  as  Hsted  below  were 
taken  from  Rel.  2  as  typical  for  double -base  propellants,  and  motor 
dimensions  were  taken  from  Ref.  3. 
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Thermal 


Component 

Material 

Density 

Ib/in^ 

Specific  Heat 
Btu  /lb-“F 

Conductivity 
Btu /in-sec- *F 

Motor  Tube 

2014  A1 
Alloy 

0.  10 

0.  22 

20.  6  X  lO"*^ 

Inhibitor 

Ethylcellu- 

lose 

0.  04 

0.  40 

2.  60  X  10"^ 

Propellant 

N-4  • 

0.  058 

0.  45 

2.  78  X  10"^ 

In  determining  temperature  histories  in  this  type  of  problem,  it  is 
the  usual  practice  to  apply  the  computed  heat  transfer  (film)  coefficients 
and  the  external  driving  temperatures  as  the  forcing  conditions  which 
cause  temperature  build-up  in  the  motor.  In  this  case,  however,  the 
measured  temperature  of  the  motor  tube  was  used  as  the  forcing  function 
on  the  grain,  thereby  isolating  the  motor  from  its  environment  for 
which  characteristics  could  not  be  accurately  determined  in  the  actual 
tests. 


ANALOG  COMPUTER  RESULTS 

Figure  1  shows  the  actual  temperature  hi.stories  of  the  motor  wall, 
the  inhibitor-grain  interface,  and  the  grain  interior  near  the  perforation 
in  curves  1,  2,  and  3,  respectively.  The  measurements  were  made  at 
identically  located  points  on  the  *op  and  one  side  of  the  motor;  however, 
only  the  grain  interface  temperatures  were  strongly  position-dependent, 
accounting  for  the  different  curves  2  and  2'.  Computed  temperature 
histories  at  the  same  positions,  assuming  grain  properties  as  given  in 
Ref.  2,  are  shown  for  zero  air-gap  by  curves  2A  and  3A  and  for  constant 
0.  045-inch  air-gap  at  the  motor  side  by  curves  2B  and  3B.  The  latter 
case  corresponds  to  the  maximum  allov/able  air-gap  under  HPAG  motor 
specifications.  These  computed  curves  for  extreme  conditions  show 
poor  agreement  with  the  measured  temperatures,  and  additional  computa 
tions  (not  included  here)  disclosed  that  neither  a  reduced  constant  air- 
gap  nor  different  grain  thermal  properties  appreciably  improved  the 
agreement. 

Figure  2,  curves  2C  and  3C,  compares  the  actual  results  with  the 
computed  temperatures  with  properties  as  before,  except  that  a  0.  020- 
inch  initial  air-gap  at  the  motor  side  continuously  changed  until  full 
contact  between  grain  and  motor  tubes  was  obtained.  The  computed 
values  for  a  grain  with  increased  thermal  di'ffusivity  but  the  same  air- 
gap  conditions  is  shown  by  curves  2D  and  3D.  A  considerably  improved 
fit  with  actual  data  is  evident  when  compared  with  Fig.  1;  however, 
reduction  of  the  initial  air-gap  and  a  smaller  increase  in  thermal  dif- 
fusivity  as  shown  in  curves  2E  and  3E  of  Fig.  3  are  required  for  an 
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Th^  ell^ct  ol^rain  heating  on  grain  growth  is  shown  for  a  typical 
case  in  Fig. '4,  andrthe  ^basis  for  computing  growth  is  described  in 
"Method  of*Calculating  Thermal  Expansion"  (p.lO). 

HEAT  TRANSITER  COEFFICIENTS 

♦;^The  over-all  heat  transfer  coefficient  for  the  heating  run,  as  cal¬ 
culated  in  Ref.  2,  is  shown  in  Fig.  5.  In  addition,  Fig.  5  shows  the 
continuous  heat  transfer  coefficients  calculated  in  this  study  by  two 
diffe'rent  methods. 

One  method  required  a  plot  of  dTj/dt  versus  time  for  the  motor 
tube  (to  avoid  differentiating  on  the  computer)  which  was  used  in  solving 
the  equation 


^1-2  +Kl 


/dtj 


T«,  -  T 


1 


where 


h  =  film  coefficient,  Btu  /in^-sec-®F 

^11-2 

C, 


heat  transfer  rate  wall-to-grain,  Btu  /in  -sec 

2 


1 


heat  capacity  of  motor  wall,  Btu  /in  -®F 
motor  wall  temperature,  *F 
oven  temperature,  °F 


The  film  coefficient  (h)  is  obtained  by  dividing  the  instantaneous  heat 
exchange  at  the  motor  wall  by  the  difference  between  the  motor-  wall 
and  the  oven  temperatures. 

The  second  method  required  manual  adjustment  of  the  film  coef¬ 
ficient  by  a  trial-and-error  operation  so  that  continuous  agreement 
between  actual  oven  temperature  and  motor-wall  temperature  was  main¬ 
tained.  Both  methods  indicate  lower  film  coefficients  during  the  early 
part  of  the  heating  period  and  support  the  belief  of  the  Ref.  1  authors 
that  oven  variables  did  cause  variations  in  the  film  coefficient. 
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FIG.  2.  Heating  Curves  for  the  5.0-Inch  HPA 
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Z.  Heating  Curves  for  the  5.0-Inch  HPAG  Rocket  Motor. 


TEMPERATURE 


FIG.  3.  Heating  Curves  for  the  5.0-Inch  HPj 
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fc.  3.  Heating  Curves  for  the  5.0-Inch  HPAG  Rocket  Motor. 
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FIG.  4.  Computed  Grain  Growth  During  Heating  Period. 


FIG.  5.  Computed  Heat  Transfer  Coefficients 
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METHOD  OF  CALCULATING  THERMAL  EXPANSION 

OF  HPAG  MOTOR  GRAIN 

Basic  Assumptions 

1.  Coefficients  of  thermal  expansion,  as  taken  from  Ref.  2,  are 
taken  equal  for  radial  and  tangential  directions. 

2.  Grain  perforation  is  taken  as  a  circular  cylinder  with  all  growth 
non-elastic  and  non-plastic. 

3.  Nominal  growth  of  each  zone  in  grain  is  proportional  to  tem¬ 
perature*  increase,  and  growths  are  additive  radially. 

4.  All  grain  deformations  are  plastic  and  are  additive  radially. 

Calcvilattons 


In  the  analog  computation  for  this  problem,  the  grain  model  was 
set  up  with  seven  concentric  zones.  The  following  discussion  is  based 
on  three  zones  in  the  interest  of  simplicity. 


The  radius  of  zone  (k)  boundaries  is  and  r  ,  and  the  me*an 

radius  of  zone  (k)  is 


r 


mk 


r 

n 


+  r 


n+ 1 
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.  It  is  now  assumed  that  the  expansion  of  the  perforation  and  zone  (c) 
does  not  appreciably  affect  computation  accuracy  for  zone  (b)  using  the 
same  method  as  for  zone  (c);  that  is, 

r'  .  =  r  .  +Ar.  +A<f'  +  iArf' 

mb  mb  4  c  Z  b 

with  an  error  of  lesd  than  2  percent  at  the  end  of  the  heating  period. 
Therefore,  the  radial  growth  of  zone  (b)  is 

.  =  2«J,<?AT  . 

b  b.  mb 

and  similarly  for  zone  (a), 

•  • 

A6 '  =  2<f  a  At 

a  a  ma 


The  maxim  vim  radius,  r'j,  of  the  grain  at  time.T  ,  is  accordingly: 


l(t)  ^4 


(1  +dAT^)  + 


=  r,  +  r.dAT.  +  za{J‘  AT  AT  ,+J  AT.  ), 

14  4  a  ma  b  mb  c  me 
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SECTION  II.  DETERMINATION  pF  WALL  TEMPERATURE¬ 
TIME  FUNCTIONS  USING  AN  ANALOG  COMPUTER 

. 

William  H.  Smith 


EDITORIAL  PREFACE 

Work  done  subsequent  to  the  preparation  of  this  paper  has  shown 
that  one  of  the  simplifying  as  sumptions  Y-that  of  one -dimensional  flow-- 
may  lead  to  significant  inaccuracies  in  the  cases  chosen  as  examples. 
Because  the  heat  flux  changes  appreciably  with  position  over  a  nose  cone, 
lateral  conduction  may  not,  in  general,  be  neglected  with  impunity. 

The  methods  described  are  reasonably  valid,  however,  for  cases  where 
the  heat  flux  is  constant  over  a  significant  area. 

GENERAL  DISCUSSION  • 

Solution  of  the  unsteady- state  heat  balance  equations  for  determin¬ 
ing  nose  cone  stagnation  temperatures  may  be  accomplished  using  an 
analog  computer.  This  section  discusses  an  analog  .computer  method 
which  was  used  for  this  purpose  and  as  an  illustrative  example. 

Prior  to  computer  evaluation  of  the  stagnation  temperature  of  a 
surface,  certain  functions  must  be  obtained.  The  shape,  dimensions, 
and  material  of  the  object  of  interest  should  be  known  since  density, 
thermal  conductivity,  specific  heat,  thickness,  and  possibly  other  param 
eters  are  used  as  constants  in  analog  computer  solutions.  The  recovery 
temperature  and  heat  transfer  coefficient  must  be  determined  as  a 
function  of  time.  Graphs  of  these  functions  versus  time  are  necessary 
in  order  to  insert  these  transients  into  the  function  generators  of  an 
analog  computer.  Methods  for  calculating  recovery  temperatures  and 
heat  transfer  coefficients  are  given  in  Ref.  4  and  5  and  in  an  informal 
technical  note.  1  Programs  have  been  prepared  for  computing  these 
functions  using  the  IBM  610  and  704  digital  conputers. 

Figures  6  an,d  7  show  heat  transfer  coefficient  and  recovery  tem¬ 
perature  as  a  function  of  tirpe  for  the  flight  profile  chosen  for  the  ex¬ 
ample.  Figures  8  and  9  show  the  two  nose -cone  configurations  which 
were  evaluated. 


U.  S.  Naval  Ordnance  Test  Station.  An  Estimate  of  the  Recovery 
Temperature  and  Heat  Transfer  Coefficient  over  the  Sidewinder  1C 
Nose  at  Mach  numbers  1.9,  2.  6,  and  3.  0,  by  L.  L.  Doig.  China  Lake, 
Calif.  ,.NOTS,  27  January  1958.  (NOTS  Technical  Note  4061-5), 
CONFIDENTIAL. 
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Certain  simplifying  assvimptionB  were  made  in  obtaining  the  solution. 
These  include: 

•  1.  One-dimensional  heat  flow 

2.  Radiation  and  convective  heat  transfer  from  the  skin  to  the 
interior  of  the  missile  are  neglected 

3.  No  ablation.. 

The  term  "REAC"  on  the  graphs  refers  to  curves  obtained  from  the 
analog  computer  solutions.  The  circuits  shown  in  Fig.  10  and  11  are 
representative  and  do  not  include  such  things  as  intermediate  amplifiers 
and  scaling  factor  potentiometers  which  are  needed  in  actual  analog 
computer  solutions.  A  thorough  knowledge  of  analog  computer  circuits 
is  needed  for  the  solution  of  particvilar  problems  and  a  person  with  such 
a  knowledge  should  be  consulted  before  any  solutions  are  attempted. 


THERMALLY  THIN  SECTIONS 

Sections  which  have  a  temperature  gradient  of  10%  or  less  were 
considered  to  be  thermally  thin  in  this  application.  The  percent  tem¬ 
perature  gradient  across  particular  designs  may  be  obtained  by  using 
Section  IV  of  this  report.  The  thermally  thin  design  considered  is  shown 
in  Fig.  8. 


The  equation  used  to  determine  the  temperature  for  thermally  thin 
sections  is 


'dT 

“dt 


w 


h 

w<i.T 


-  T  ) 
w 


(1) 


The  computer  circuit  (Fig.  10)  needed  to  solve  this  equation  is  quite 
simple.  Approximately  4  hours  were  needed  to  wire  the 'control  panel, 
put  the  recovery  temperature  and  heat  transfer  coefficient  into  the  function 
generators  of  the  Mid-Century  Analog  Computer,  and  read  out  solutions 
for  the  problem. 

Figure  11  shows  the  skin  .temperature  as  a  function  of  time  using  re¬ 
sults  obtained  from  computer  solutions  of  Eq.  1.  The  values  of  emis- 
sivity  were  arbitrarily  chosen  at  0.  0,  0.  5,  and  0.  85. 

Several  methods  of  numerical  integration  are  available  for  hand 
solution  of  Eq.  1.’  These  include  repetitive  calculations,  trial  and 
error,  and  the  Runge-Kutta  methods,  and  are  fully  explained  in  Ref.  6. 
When  the  radiation  effect  is  included,  the  only  one  of  these  methods  ap¬ 
plicable  is  the  trial-and-error  solution.  However,  the  effects  of 
radiation  may  normally  be  neglected  below  1000°F. 
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FIG.  11.  Wall  Temperature  Versus  Time--Stainless  Steel 
6.9-inch  Diameter  Dome,  0.030-inch  Thick. 

THERMALLY  THICK  SECTIONS 

If  the  temperature  gradient  across  &  particular  configuration  is 
greater  than  10%,  the  configuration  is  termed  thermally  thick.  In 
thermally  thick* designs,  conduction  must  be  considered  to  obtain  correct 
solutions.  The  thermally  thick  design  for  this  example  was  an  Astrolite 
hemisphericaf  dome  as  shown  in  Fig.  9-  As  with  thermally  thin  config¬ 
urations,  radiation  effects  on  the  outside  surface  temperature  may  be 
neglected  when  the  temperature  of  this  surface  is  below  1000  °F. 

The  general  equation  for  unsteady- state  conduction  as  given  by 
Fourier's  general  law  of  heat  conduction  can  be  written: 

(2) 


dt 


K 


2  2  2 

a'^T  ^  a'^T  ^  a'^T 

— -r  + 

5  X  ay  a  z 
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For  the  one -dimensional  heat  flow  which  is  being  considered,  the 
equation  reduces  to:  * 

^  ^  r ^'^1  /■»» 

jr  - 

Partial  differential  equations  such  as  Eq.  2  and  3  are  not  directly  sol¬ 
uble  using  analog  computers.  The  solution  for  wall  temperatures  of 
thermally  thick  sections  may  be  closely  approximated  using  the  following 
method. 


Assume  that  instead  of  a  section  with  a  continuous  temperature 
gradient,  a  model  is  devised  composed  of  a  numbqr  of  thin  slabs.  Each 
slab  has  a  uniform  temperature  across  its  width.  To  maintain  the  same 
total  gradient  from  outside  to  inside  surface  as  the  original  section  (with 
continuous  gradient),  the  temperature  gradient  of  the  model  must  take 
the  form  of  a  step  function,  each  step  taking  place  at  a  slab  interface. 
Using  such  a  model,  a  heat  balance  equation  may  be  written  for  each 
slab.  The  equation  of  the  outside  surface  slab  of  such  a  model  is: 


h 

wCdx 


(T  .  - 
r 


T  ) 
w 


wCAx 


K 


wC(Ax) 


(T 


w 


-Ti)  = 


AT 


M. 


At 


(4) 


The  heat  balance  equation  for  successive  slabs  is: 
dT, 


1 


K 


dt  ~  —r~iA  7^  ”  ‘‘*1  *2' 


dT2 

-W 


wC(Ax) 
K 


(T..  -  2T,  -  T,) 


wC(Ax) 


2  <•'^1  -  2^2  - 


(5) 

(6) 


Equations  such  as  Eq.  5  and  6  can  be  written  for  each  slab  until 
the  inside  surface  is  reached.  Since  this  method  is  at  best  an  approxi¬ 
mation,  the  effects  of  radiation  and  convection  are  usually  neglected  due' 
to  their  small  values. 

The  number  of  slabs  into  which  a  particular  configuration  is  divided 
depends  on  the  accuracy  desired  and  the  size  of  the  computer  available. 
The  maximum  permissible  temperature  gradient  per  slab  thickness  was 
set  at  10%.  Smaller  maximum  permissible  temperature  gradient  might 
be  chosen  if  greater  accuracy  is  desired. 

The  analog  computer  circuitry  needed  for  computer  solution  of 
Eq.  4,  5,  and  6  is  shown  in  Fig.  12.  This  circuit  includes  solutions 
for  the  outside  surface  and  the  first  two  slabs.  Additional  circuitry 
may  be  added  to  this  basic  circuit  until  a  model  with  the  correct  number 
of  slabs  is  obtained. 
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A  method  has  been  worked  out  for  obtaining  accurate  solutions  of 
Eq.  4  using  hand  calculations  (Ref.  7),  As  with  the  thermally  thin  case, 
however,  the  method  is  time  consuming  and  tedious  for  most  applica¬ 
tions.  Although  no  provision  for  radiation  effects  is  included  in  this 
method,  an  approximation  of  its  effect  on  outside  surface  temperature 
may  be  made  by  calculating  the  heat  lost  due  to  radiation. 

Figures  13,  14,  and  15  show  some  of  the  solutions  obtained.  These 
solutions  include  the  effect  of  radiation  on  outside  surface  temperature 
(Fig.  13),  temperature  gradient  across  the  section  transient  (Fig.  1^), 
and  the  effect  of  total  thickness  on  the  inside ‘surface  temperature  (Fig.  15). 


FIG.  13.  Outside  Wall  Temperature  Versus 
Time  for  Varied  Emissivity  (€ )- - Astrolite . 
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FIG.  14.  Wall  Temperature  Versus  Time;  Out'side,  Middle,  Inside 
Astrolite  0.  250-inch  Thick  Wall,  £  =  0.  85 


FIG.  15.  The  Effect  of  Wall  Thickness  on  Inside  Wall  Temperature 
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COMPARISON 

As  a  check  on  the  accuracy  of  the  computer  solutions,  the  method 
of  Ref.  7  was  also  used  to  obtain  solutions.  These  calculations  were 
carried  out  for  time  increments  of  1,  2,  and  3  seconds.  Figure  16  shows 
\he  results  of  these  calculations  and  the  computer  solution  at  the  lower 
recovery  temperatures  of  the  flight  profile  assumed.  The  agreement 
between  the  various  solutions  is  very  good.  Figure  17  shows  the  results 
obtained  at  the  higher  recovery  temperatures  with  radiation  not  in¬ 
cluded.  The  computer  solution  lies  generally  between  the  2-  and  3-sec 
increment  solutions. 


FIG.  16.  Temperature  Versus  Time  Comparison. 
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» 


TIME,  SECONDS 


FIG.  17.  Outside  Wall  Temperature  Versus  Time  Comparison-- 

Zero  Emissivity,  Astrolite. 

Figure  18  is  a  comparison  of  computer  and  hand  solutions  for  out¬ 
side  surface  temperature.  The  farter  response  of  the  computer  solution 
is  effectively  shown  on  this  graph. 
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FIG.  18.  Hot  Wall  Temperature  Versus  Time--Astrolite, 
6.9-inch  Diameter  Dome,  1/4-inch  Thick; 
Emissivity  =  0.  85. 


SUMMARY 

The  use  of  an  analog  computer  to  obtain  solutions  of  the  outside 
surface  temperature  of  aerodynamically  heated  bodies  results  in  large 
time  savings  over  hand  solutions  of  the  same  problem.  After  the  control 
panel  has  been  wired,  and  the  recovery  temperature  and  heat  transfer 
coefficient  functions  have  been  put  into  the  computer,  solutions  for  a 
large  number  of  emissivities ,  wall  thicknesses,  and  even  different 
materials  can  be  obtained  in  a  very  short  time.  These  changes  are 
normally  accomplished  by  changing  one  or  more  potentiometers  in  the 
computer.  A  completely  new  set  of  hand  calculations  is  required,  how¬ 
ever,  for  each  new  set  of  initial  parameters.  Also,  while  radiation 
effects  are  included  quite  simply  in  computer  solutions,  they  add  very 
much  to  the  complexity  of  hand  solutions. 

The  accuracy  of  the  computer  solution  appears  to  be  very  good, 
although  further  experimental  work  is  needed  to  show  exactly  how  good. 
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NOMENCLATURE 


Btu 


C  Specific  heat  of  material, 

h  Heat  transfer  coefficient,  — 5 — - 

*  in  -sec-*F 


K  Thermal  conductivity ,t- 


Btu 


in-sec- 

T  Absolute  temperature,  ®R 
t  Time,  seconds 
w  ‘Specific  weight,  Ib/in^ 

A  X  Slab  thickness,  inches 
Emis  sivity 


e 

cr 


Stephan  Boltzman  constant,  — « — r 

in  -sec- "R 


T  Total  thickness  in  thermally  thin  configurations, 


Subscripts 


w  Wall  conditions 
r  Recovery  conditions 
n  Slab  conditions 


inches 
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SECTION  HI.  WHEN  IS  A  SKIN  THERMALLY  THIN? 

By 

Jerome  R’.  Katz 


INTRODUCTION 

« 


In  the  problem  of  aerodynamic  heating  it  has  become  desirable  to 
know  when  a  layer  or  skin  of  an  object  flying  in  air,  and  hence,  sub¬ 
jected  to  heating,  can  be  considered  thin.  That  is,  when  can  the  tem¬ 
perature  be  considered  uniform  throughout  the  skin?  Thermodynamically, 
a  skin  can  be  considered  thin  when  the  material  has  a  high  conductivity, 
or  when  the  time  for  which  it  is  subjected  to  heating  is  long.  However, 
for  this  report  a  thin  skin  will  be  defined  in  terms  of  the  thickness  of 
the  material,  which  will  be  a  single  layer  metal.  The  skin  is  approx¬ 
imately  thin  when  the  derivative  of  temperature  with  respect  to  the 


thickness  can  be  considered  zero  i 

dx 


temperature  gradients 
A  method  for  determining  a 


across  its  thickness  can  be  neglected, 
skin  thickness  is  discussed  in  this  section  for  1-,  5-,  and  10-percent 
variation  in  temperature  throughout  the  material.  The  calculations  for 
heat  transfer  were  performed  for  a  flat  plate  with  heat  being  added  on 
one  side  and  insulation  on  the  other  (Ret.  7).  For  the  purpose  of  this 
report,  the  insulation  will  simulate  a  no-heat  transfer  condition. 


METHOD 

* 

Reference  7  contains  values  of  transient  temperatures  which  have 
been  normalized.  The  normalized  temperatures  are  represented  by 
u  and  are  solutions  of  Eq.  1  of  Ref.  7  tabulated  for  a  given  Biot  modulus, 
OL ,  x/L  ratio,  and  Fourier  modulus,  Z  .  The  Fourier  modxilus,  T  ,  u,  and 
the  Biot  modulus,  ct  ,  are  defined  as  follows: 


cx 


hL 

K 


U  =  (U  -U.)u+U.  ,  (0:Su^l) 

O  1  1 

kL^T  =  t,  (O.'OOl^  T  £  1, 000) 

Since  one  side  was  subjected  to  heat  transfer,  Eq.  1  can  be  written: 

h  L 

■  ^ 

o  TT 

where  L  repres&nts  plate  thickness.  The  equation  representing  the 
conditions  on  the  other  side  (no  heat  transfer)  is  defined  by  a..  =  0. 

• 

28  .  •  • 


NAVWEPS  REPORT  7653 


Each  page  of  the  tables  ih  Ref.  7  contains  the  above  data  for  a  different 
value  oi^o  ranging  frona  0.  001  to  1, 000  thusly:  =  0.  001,  0.  002, 

0.004,  0.006,  0.01,  0.02,  0.04,  0.06,  0.1,  etc.  °Hence,  this  gives 
25  pages  of  data.  Now  for  each  page,  or  each  value  of  <ac  ,  a  value  of 
T  can  be  foiind  corresponding  to  a  1-,  5-,  and  10-percenP  variation  in 
temperature  in  the  following  manner: 


First  take  the  ratio  rr  for  each  T  entry.  By  linear  interpola- 

Uo 

tion,  when  the  value  of  the  above  ratio  was  found  to  be  0.  01,  0.  05,  and 
0.  10,  a  similar  interpolation  process  was  used  to  yield  corresponding 
values  oft.  This  gives  for  each  at  ,  a  T  for  which  the  temperature 
variation  is  1 ,  5,  and  10  percent  throughout  jhe  skin. 


Next,  a  was  plotted  as  a  function  of  T,  =  fj  {^)  with  the  1-,  5-, 
and  10-percent  temperature  variation  as  parameters.  Now  we  make  use 
of  Eq.  3  and  4  of  Ref.  7  and  the  following  relation: 


Combining  these  expressions  yields  another  function  of  in  terms  of 
T  where  L  has  been  eliminated. 


Let  C  = 


c^K 


o 

Then  Eq.  5  (Ref.  7)  can  be  written  as: 

a  2  _  C 
o  T 


For  a  given  C,  a  is  a  function  of  T.  Values  of  C  ranging  from  0.  0001 
to  10,  000  and  T  from  1  to  50  were  chosen  in  order  to  include  a  wide 
range  of  physical  conditions.  Then  «  was  calculated  as  a  function  of 
Z,  a  =  f2  (T)  on  the  IBM  704.  This  °data  along  with  the  previous 

(Z)  curve  were  plotted  in  Fig.  20  for  the  range  of  C  values.  The 
points  of  intersection  of  the  two  functions  yield  values  of  iZq  that  when 
substituted  into  the  equation 


a 

o 


h  L 
o 


-  ~K- 


will  yield  values  of  L.  for  which  the  skin  of  an  object  of  given  material 
can  be  called  thin  within  the  range  of  values  of  C  and  a  given  percentage 
of  temperature  variation. 


The  author  gratefully  acknowledges  the  assistance  of  Gladys 
Radeck  for  IBM  704  calculations. 
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FIG.  19.  Biot  Modulus  as  a  Function  of  T  for  1-,  6-,  and  10 
Percent  Temperature  Variation  in  a  Flat  Plate. 
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EXAMPLE 

To  illustrate  the  method  described  in  the  text,  a  hypothetical 
problem  is  selected.  It  is  desired  to  know  for  a  .1%  temperature  vari¬ 
ation  when  a  flat  plate  of  aluminum  alloy  can  be  called  thermally  thin, 
after  being  subjected  to  heating  for  10  seconds.  First  the  physical  and 
thermal  properties  are  found. 


=  169 

Ib/ft^ 

K 

=  117 

Btu/hr-ft-'F 

c 

=  0.  18 

Btu/lb-"F 

f 

=  10 

sec 

Selecting  a  maximimi  heat  transfer  coefficient  h^  =  0.  01  Btu/ft  -sec 
(Ref.  8),  the  values  of  C  can  be  evaluated. 

c  -  **^0^ 


(10)  (0.  01)“ 


.  (3600) 


3.  6 

C  =  0.  001 

Referring  to  Fig.  20,  the  line  corresponding  to  C  =  0.  001  is  used, 
and  the  intersection  of  this  curve  with  the  1%  temperature  variation 
curve  is  found.  This  yields  an  -  0.  0047.  Now  using  Eq.  1  (Ref.  7), 
the  desired  thickness  can  be  obtained. 


CL  = 

o 


(0.  0047) (0.  117) 


0.  184  inches 


Approximately  the  same  thickness  can  be  obtained  by  using  an  average 

heat  transfer  coefficient  h  of  0.  003  Btu/ft^-sec- “F  and  a  time  of  100  se 

o 
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NOMENCLATURE 
thp^ 


cp 

c  Specific  heat  of  object,  Btu/lb-*F 

2 

h  Heat  transfer  coefficient,  Btu/ft  -sec-*F 

K  Thermal  conductivity,  Btu/hr-ft-*F 

ij  _  f  =  IT  = 

thermal  diffusivity  —  ~  "K” 

^  c/> 

L  Total  thickness  of  skin  inches 

t  Time,  sec 

U  Temperature,  *F 

u  Normalized  temperature  from  tables  of  Ref.  7 

X  Distance  through  the  skin,  inches  from  side  exposed  to  heating 

a.  Biot  modulus  based  on  skin  thickness 

.3 


Density,  Ib/ff 
T  Fourier  modulus. 


t 


kL 


Subscripts 

o  Conditions  at  the  exposed  side 
i  Conditions  at  the  insulated  side 
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